Purification to Locally Maximally Entangleable States 
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Locally maximally entangleable states (LMESs) constitute a large set of multipartite states, con- 
taining for instance all stabilizer states. LMESs are uniquely characterized by (2™ — 1) phases, where 
n denotes the number of qubits. We consider here those LMES whose phases are either or 7r and 
present a multipartite entanglement purification protocol for arbitrary such states. In contrast to all 
previously known recurrence protocols this protocol uses a novel ingredient, which is required due to 
the quantum correlations contained in the various LMESs. We compare this scheme to previously 
known entanglement purification protocols and show that the direct purification performs better 
than previously known protocols. 

PACS numbers: 



I. INTRODUCTION 

The focus on quantum entanglement has shifted from a 
puzzling phenomenon of fundamental interest to a valu- 
able resource in the context of quantum information pro- 
cessing. Several applications of entanglement are known, 
ranging from quantum metrology over quantum commu- 
nication and computation to (multiparty) security appli- 
cations, see e.g. [1-4]. In all these applications pure 
states, often distributed among several spatially sepa- 
rated parties, are required. Despite spectacular exper- 
imental progress, the creation and maintenance of mul- 
tiparticle entangled states remains a difficult task. Noise 
and decoherence limit the local generation of states, 
while channel noise is a main obstacle when creating dis- 
tributed entanglement. 

Entanglement purification is known to offer a possi- 
ble way to overcome these limitations [5] . Entanglement 
purification protocols were introduced for bipartite states 
[5], where it was shown that despite significant noise, e.g. 
resulting from sending particles through noisy channels, 
several copies of noisy entangled states allow one to gen- 
erate fewer copies with high fidelity. Later entanglement 
purification protocols were introduced and analyzed for 
certain multipartite entangled states, in particular graph 
states or stabilizer states [6-12] and W-states [13]. The 
possibility to purify certain states with the help of bipar- 
tite or multipartite entanglement purification protocols 
ensures that these states can be generated or maintained 
with high fidelity, even in the presence of noise. This is 
important when using such states as resources for quan- 
tum information processing or other tasks. In addition, 
its also a question of fundamental interest for which types 
of quantum correlations a purification process is possible. 

Moreover, studying certain classes of states and inves- 
tigating their entanglement features and possible appli- 
cations has proven a to be a successful approach towards 
the understanding of multipartite entanglement -a com- 
plex problem that is still far from being completely un- 
derstood. Graph states or stabilizer states can be viewed 
as such an example. They constitute a large class of 
multiparticle entangled states, including several interest- 



ing and highly entangled quantum states such as the 
Greenberger-Horne-Zeilinger state, codewords for error 
correcting codes or the 2D cluster state [14]. The lat- 
ter serves as a resource for measurement based quantum 
computation [4, 15]. 

Here we introduce a novel entanglement purification 
protocol for an even larger class of multiparticle entan- 
gled states, the the so-called locally maximally entan- 
gleable states (LMESs) [16], thereby paving a way for 
possible practical applications of such states in the pres- 
ence of channel noise and decoherence. LMESs are de- 
fined as those states where all quantum information can 
be maximally washed out by coupling each particle to a 
local auxiliary qubit and are equivalent to those states 
that can be used for optimal local encoding of classi- 
cal information. LMESs correspond (up to local basis 
change) to a coherent superposition of all basis states, 
where all coefficients have equal norm and therefore are 
characterized by a phase. This leads to a total of (2™ — 1) 
parameters describing such a n-qubit state. In this pa- 
per we concentrate on LMESs where all phases are either 
or 7r, i.e. all coefficients are cither (+1) or (-1). In 
this sense LMESs can be viewed as a generalization of 
graph states [17, 18] or stabilizer states, where however 
stabilizers are no longer simple tensor products of Pauli 
operators. For all such LMESs, we provide a so-called re- 
currence protocol that allows one to purify a mixed state 
with sufficiently high initial fidelity to any desired LMES. 
More precisely, we present protocols that operate on two 
identical copies and which are capable of increasing the fi- 
delity provided the initial fidelity is sufficiently high. The 
purification process is divided into several sub-protocols. 
Each sub-protocol provides purification only with respect 
to certain parties, and the different sub-protocols need to 
be combined to achieve an overall entanglement purifica- 
tion with respect to the desired state. 

The approach is similar to the one presented for the 
purification of graph states [8-11], and in fact the result- 
ing purification maps with respect to diagonal coefficients 
in the LMES basis turn out to be equivalent. However, 
new ingredients not used so far in entanglement purifica- 
tion protocols are needed to achieve this goal. It turns 
out that controlled-not operations or local stabilizer mea- 
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surements [12] are no longer sufficient, and new types of 
parity measurements are required. In contrast to Bell 
state or graph states purification protocols, it is also not 
possible to bring initial states to some standard form via 
local depolarization, i.e. diagonal in a LMES basis (which 
we call LME diagonal) as we show in this article. Nev- 
ertheless, we demonstrate that entanglement purification 
is possible for a large class of noisy LMES. 

This paper is organized as follows. In Sec. II we intro- 
duce our notation and review graph states and LMESs. 
In Sec. Ill we introduce and analyze our entanglement 
purification protocol. We first focus on a specific class of 
LMESs, and then show how this approach can be used to 
purify arbitrary LMESs with phases and it. In Sec. IV 
we discuss depolarization to LME diagonal states. There 
we show the impossibility of local depolarization, and 
provide an example where depolarization may spoil pu- 
rification. In Sec. V we analyze several examples and 
compare our direct multiparticle entanglement purifica- 
tion protocol with possible alternative approaches, where 
we show an advantage of our approach in certain param- 
eter regimes. We summarize and conclude in Sec. VI. 

II. STABILIZER STATES, LOCALLY 
MAXIMALLY ENTANGLEABLE STATES 

In this section we introduce our notation and review 
some basics concerning stabilizer and locally maximally 
entangleable states (LMESs). 

A. Notation 

Throughout the manuscript we use the following nota- 
tion. X, Y, Z denote the Pauli operators and 1 denotes 
the unnormalized identity operator. The subscript of an 
operator (state) denotes the qubit the operator (state) is 
acting on (describing) resp., e.g. X\ denotes the X oper- 
ator acting on qubit 1 and pi denotes the reduced state of 
qubit 1. Unless stated otherwise, a superscript of an op- 
erator denotes the exponent of the operator, e.g. X° = 1, 
and X 1 = X. Moreover, we will often us multi-indices, 
which we denote by bold letters, like k = (fci, . . . , k n ). 

B. Stabilizer states and Graph states 

A n-qubit stabilizer state, |^), is uniquely defined via 
a set of n commuting operators {Si, S„} in the Pauli 
group [3]. The state |\&) is the unique eigenstate to eigen- 
value + 1 for all the operators Si, i.e. Si |<£») = |<&) Vi iff 
|$) = |*). The group generated by {Si, S n } is called 
stabilizer group of We will call the generators, Si, 
the stabilizer of I*). Since the eigenvalues of each of the 
Pauli operators are ±1 the eigenvalues of any stabilizer 
are also ±1 and are both 2™ _1 — fold degenerate. Sim- 
ple examples for stabilizer states are Bell states, e.g. the 



stabilizer of |$+) = 1/V2(|00) + |11)) is generated by 
{Si =X®X,S 2 = Z®Z}. 

Note that each stabilizer state is up to local Clifford op- 
erations equivalent to a graph state [17, 18] . Graph states 
can be associated to a mathematical graph G = (V,E) 
consisting of a set of vertices V — {1, . . . , n} and edges, 
E C € V, i 7^ j}. The graph state corre- 

sponding to a graph (V, E) is = U p h |+)® n with 
U P h = Udj^EUij. Here, U l0 = |0> (01,® + (l^Zj 
denotes the 2— qubit n— phase gate. The neighborhood of 
qubit k, Nk := {/ € V\{k, I) € E}, is the set of all vertices 
which arc connected by an edge to qubit k. It is easy to 
see that the stabilizer of a n— qubit graph state is gener- 
ated by S t = UphXiU^ = X, <S) jeN . Zj, for i e V. Note 
that the common eigenbasis of the stabilizer operators, 
Si, is given by {|* k ) = Z fcl ®Z k *® . . .®Z k - |*)} fcie{0 ,i}, 
with l^o) = I*)- Obviously, all elements of this ba- 
sis are stabilizer states and it is easy to verify that 
<%|*k) = (-l) fei |*k)- Thus, the index k contains all 
the information about the eigenvalues of the stabilizers 
to which l^k) belongs to. 

A useful concept in the context of graph states is col- 
orability. Qubits which are not connected among each 
other by a phase gate can be grouped into a set called 
color. A graph state is said to be k— colorable if one 
can group the set of vertices into k subsets (colors) such 
that there are no edges among qubits within one sub- 
set. For instance, the 3-qubit GHZ state, \GHZ) = 
l/v / 2(|000) + |lll)) is up to local unitary operators equiv- 
alent to U12U13 |+)® 3 and is therefore 2-colorable. Note 
that local unitary equivalent graph states might have a 
different colorability and that the minimal number of re- 
quired colors is in general not known. 

As mentioned in the introduction multipartite entan- 
glement purification protocols have been derived to pu- 
rify to arbitrary graph states [7-10]. Here, we briefly 
summarize this purification protocol. Let us consider 
as a target state a 2-colorable n-qubit graph state, 
|*). The aim is to transform N copies of an input 
state p into via LOCC. Since the set {|5'k)}ke{o,i}" 
forms an orthonormal basis, p can be written as p — 
Sk l l*k) (^ll- It has been shown that any input 
state can be transformed via LOCC into a state diagonal 
in a graph-state basis, p out — X)k ^ k l* k ) (*k|, without 
changing the fidelity F = (&\p\V) [24]. Since all states 
|\l/k) are local unitary (LU) equivalent to each other, the 
aim of entanglement purification is to gain information 
about the index k, i.e. the eigenvalues of the stabiliz- 
ers. Similarly to the recurrence protocol for Bell state 
[5], the main idea is to transfer some information about 
the eigenvalues of the stabilizers from the first copy of the 
input state, pi to the second copy, p2- This is achieved 
via multilateral CNOT operations. A measurement on 
the second copy reveals this information, which allows to 
purify the remaining copy. Since the qubits belonging 
to the same color do not interact, the information about 
the eigenvalues of the corresponding stabilizer can be ob- 
tained simultaneously. Later on, it has also been shown 
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how arbitrary fc-colorable graph states can be purified by 
using 2-colorable graph states as auxiliary states. The 
required 2-colorable graph states can be obtained from 
the fc-colorable graph states via LOCC. As we will see, 
the tools required for the purification of LMESs will be 
significantly different. In particular, we will derive an en- 
tanglement purification protocol where other operations 
than the CNOT gates are required to reveal the desired 
information. 



C. Locally Maximally Entangleable states 

In this section we review some properties of the class 
of LME states for which we derive purification protocols 
[16]. In the following we will consider n-qubit states and 
denote by V = {1, . . . , n} the set of all qubits. As men- 
tioned in the introduction, a state is called LMES if it is 
possible to attach (in a certain way) to each system qubit 
a local auxiliary qubit such that the whole information 
of the system is washed out. That is the reduced state of 
the n-system qubits is completely mixed after attaching 
the auxiliary qubits. It has been shown that a n-qubit 
state is LME iff it is local unitary (LU) equivalent to a 
state of the form 

\*) = u 1 ,..., n l[u iki ,... !iki ...l[u i \+f n , (1) 

where U ikii ... tikm = l+(e^-l) |1) (l\® m denotes a phase 
gate which is acting on the m qubits, {ik i: ■ ■ • ,ik m } re- 
spectively. We call such a phase gate a pure phase gate 
of order m. Equation (1) reveals a way to generate the 
state, namely by applying the appropriate phase gates to 
a product state |+)® n . Looking at this generation pro- 
cess one can see that graph states are a subclass of LME 
states. They are exactly those states which arise if one 
uses only 2-qubit interactions with n phases. Notice that 
LMESs can be used to generate an arbitrary state of an 
even larger class of so-called M-states, which are used in 
the context of classical simulation of quantum computa- 
tion [19]. It is straightforward to show that any n-qubit 
M-state can be obtained from a n + 1-qubit LMES by 
performing a measurement on a single qubit. 

Due to the definition of LMESs one can generate an 
orthonormal basis of the Hilbert space by applying in- 
dependent local unitary operations to the state. In par- 
ticular, for an LMES \*5f) as given in Eq.(l) one can ap- 
ply local Z operations to generate an orthonormal basis 
{|* k ) = Z kl <g> Z k2 <g> . . . <g> Z k " |*)} fcie0 ,i. Moreover, it is 
easy to verify that the group generated by the operators 
Si = U ph XiUl h stabilizes the LMES |#j = U ph |+)® n . 
Like in the case of stabilizer states the operators fulfill 
the conditions Sf = 1, [Si, Sj] =0 and Si = Sj, for any 
i,j G V. Note however, that in contrast to stabilizer 
states the stabilizer of a LMES are no longer necessarily 
elements of the Pauli group, but act non-locally on sev- 
eral qubits. Analogously to the stabilizer states, it can be 



easily verified that Si |*k) = (-l) fei l^k)- Hence, | \& k > is 
an eigenstate of Si to eigenvalue (— l) fe \ 

Similarly to graph states we introduce the notion of 
neighborhood and colorability for LMESs. A qubit i is 
said to be in the neighborhood of qubit j if U p h con- 
tains a pure phase gate acting on both of them. The 
neighborhood, i.e. the union of all neighbors, of qubit 
j will be denoted by Nj. We group those qubits of 
a LMES which are non-interacting among each other 
in different sets A = {a\, . . . aj A }, B = {b\, . . . , bi B }, 
etc for some Ia,Ib,--- € {l,...,n}. That is, there is 
no interaction between any two qubits belonging to the 
same set, say A. In order to ease the notation we de- 
note by A, B, etc. both, the set of qubits belonging 
to that color as well as the color itself. A LMES is 
then called fc— colorable if k is the number of required 
colors. Whenever a /c-coloration of a LMES has been 
chosen it will be advantageous to divide the multi-index 
k into the k multi-indices k^ks, . . . corresponding to 
the different colors. We will then for instance write 
|*k A ,k B ,...) instead of |*k) when referring to the state 

From now on we will restrict ourselves to 7r-phases, 
i.e. any pure phase gate is of the form U — 
diag(l, 1, ... ,1,-1). Note that this choice is not only 
for the sake of mathematical simplicity but also because 
of physical reasons as this class seems to be the most nat- 
ural generalization of graph states. We call a fc-colorable 
LMES, [*£) = Uph |+) ™, regular if the decomposition of 
Uph into pure phase gates contains only pure phase gates 
of order k. Note that this implies that all pure phase 
gates are acting on one qubit per color. An example of 
a regular 3-colorable LMES is U123U234U345U456 |+)® 6 
with A = {1,4}, B = {2,5}, and C = {3,6}. Note 
that for a regular LMES, the set of neighbors of qubits 
in some color, say A, Na = Ua-€A-^<»i' contains all 
qubits which are not in A. That is Na — N\A = A. 
Thus, for any regular LMES, |*) = U ph |+)®", we have 

U P h = U a ^AuiT' Na '\ where U^ ] = |0) Oi (0| ® 
1 + |1) . (1| ® U ai denotes the phase gate (not neces- 
sarily pure) which entangles qubit with all its neigh- 
bors [25]. For instance for the example above we have 
y|2,{i,3,4}) = Ul2aU2M = | ) (0 | 2 g, i + |i) (i| 2 U 13 U 34 . 
We will also use that U p h — J2k \^)a (^1 ® Uk, where 
Uk = U k ^ ■ ■ ■ Ua^l is acting on the neighbors of A. Us- 
ing that U ai = f/J it is easy to verify that the stabilizers 

S at are given by S ai = U ph X ai U^ h = X ai <g> U a% - More- 
over, since U ai is not acting on any qubit in A for any 

i, we have that S J a \ ■ ■ ■ Sl\*\ = {X A ) iA <g> U- y In the ex- 
ample above we have for instance, Na — {2,3,5,6} and 
Uph = ]C k l k ){i 4} ( k l ® Uk, where U o = 1, U i = U 4 = 
U2 3 U 35 U 56 , U w = U x = U 23 and U n = UJJ A - U 35 U 56 
and Si=X 1 ® U 23 - 
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III. PURIFICATION PROTOCOL 



Using the same notation as before, our aim is to purify 
an arbitrary n-qubit state, p, to a regular n-qubit LMES, 
\^>). Without loss of generality we take |\&) = |^o), which 
is considered to be the element of the LME basis with 
which p has the largest overlap. In the last subsection 
of this section we address the problem of purifying to 
arbitrary 7T-LMES. 

As usual in recurrence protocols, one assumes an ar- 
bitrarily large number of noisy copies at one's disposal. 
Since the parties are spatially separated from each other, 
they are restricted to LOCC operations. 

Let Pa denote the projector onto the subspace spanned 
by those basis states, for which Ua = 0, i.e. Pa = 
Skj 1*1^=0,0 (*k A =o,kj|- First of all, we show how 
one can transform two copies of a state p into one n-qubit 
state, p with tr(P^p) > tr(PAp) in case ti (Pap) was suf- 
ficiently large. That is we purify color A. Then we will 
combine the purification protocol of all colors in order to 
increase the fidelity, F = (*| p |*) = tr(P^P B • • • P N p). 

We consider as our target state a fc-colorable regu- 
lar LMES, |*) = U ph 1+)®". Using that the states 
l^k) = Z kl ® Z k2 <g> . . . <g> Z kn |*) form an orthonor- 
mal basis, we express the input state in this basis, 
p = J2k.\ l*k> We derive now an LOCC protocol 
which achieves the following mapping: 



|*k A ,k A ) |*u, U > |*k A ,k A ©i A ) h®\ a \ \k A e \a) , (2) 



where the last \A\ qubits are auxiliary qubits, which 
will be measured in the A -basis. Only if all measurement 
outcomes arc +1 the state will be kept. Otherwise the 
copy is discarded. Thus, this protocol enables one to 
measure the parity ® 1^. Hence, in case the initial 
weight tr(PAp) was large enough, the probability of being 
in the right subspace with respect to color A, i.e. k^ = 0, 
after the first purification step, is increased. 

Considering Eq.(2) it is easy to see that an 
arbitrary input state p = ^ k k , 
will be transformed 



'k A ,ki 



y A£>k) (*k'| 

to the unnormalized state 
Ma) (*k' A ,k^ei' A ■ In particu- 
lar, for an LMES-diagonal state, p = ^ k Ak l^k) (^k| 

we obtain £ \k A ,k A \k A ,i A \^k A ,k A ®i A ) (^k A ,k A & A \- 
The whole purification protocol consists then of a 
sequence of purifications of all colors. Even though the 
purification of some color, B, introduces again some 
noise in color A, i.e. by purifying color B one might 
decrease tr (P^p), an overall increase of the fidelity, F is 
achieved. In particular, the target state, is a fixed 
point of this protocol. 

Let us now explain in detail how the purification pro- 
tocol for one color, say color A, works. 



A. The purification protocol for one color for 
regular LMESs 

Let us denote by aj, A 1 (a 2 , A 2 ) the qubits and the set 
of qubits belonging to color A of the first (second) copy 
resp. Then, the purification protocol for color A, denoted 
by Va, reads as follows. 

i) All parties holding a qubit of color A apply locally a 
CNOT-gate, U c = |0) (0|<g>l+|l) (l|<g>X, between 
qubit aj and af. In total the unitary operation 

2 1 

Ua,eA U c " a% is applied. 

ii) Each party holding a qubit in the neighbor- 
hood of A, x e Na, creates locally the max- 
imally entangled two-qubit state, l^ -1- ),,,/ x „ and 
measures qubit x x ,x 2 and x' in the GHZ- basis, 
{A <g> 1 <g> Zi \GHZ) , 1 <g> X <g> Zi \GHZ) , 1 <g> 1 <g> 
X l Zi \GHZ)}ij=QA- Only those outcomes where 
one of the states 1 12 <g> W» \GHZ) with W l <G 
{A J Z k }j t k=o,i has been measured are kept and the 
operation W* is applied to qubit x" . Otherwise the 
copy is discarded. 

As we will show below, after the first two steps the 2n- 
qubit state |*k A ,k A ) |*U,U ) wm t* e transformed into the 
(n + |A|)-qubit state \^k A ,k A ®\ A ) H®\ A \ |k A 8 U) (we 
ignore here those qubits which have been measured in 
the GHZ-basis and therefore factorize). Note that to 
each party in color A one of the last \A\ qubits belongs 
to (see Fig. 1 for an illustration). More precisely, qubit 
a 2 is after steps (i) and (ii) is in the state H \ k l A © l A ). 

iii) All parties holding a qubit of color A measure their 
qubits af in the A-basis. They keep only those 
instances where the outcome +1 is obtained and 
discard the rest. Thus, only those copies, where 
kyt ® I a = are kept. 

Before we show that this protocol achieves indeed 
the desired transformation, we want to stress that in 
contrast to existing recurrence protocols, we need here 
the projection onto the GHZ-state (step (ii)). That 
is, in contrast to the Graph state protocol, it is not 
enough to apply multilateral CNOT-gates. The rea- 
son for this is the different kind of quantum correla- 
tions contained in LMES. To make this statement more 
precise let us consider the simple example of the three- 
colorable regular LMES, |#) = U 123 |+)® 3 . Consider- 
ing the application of U^ 1 l^k)^ l*i)i'2'3' we obtain, 

^ lil ^i'23|*k)i 23 l*i)i'2'3'- Thus > thc CNOT gate en- 
tangles the two copies in such a way that no other local 
CNOT gate can disentangle them again keeping some in- 
formation about the eigenvalues of S ai in the second copy. 
Therefore, any subsequent measurement on the second 
copy, which is required to "project" the remaining copy 
onto the desired subspace, will destroy the entanglement 
contained in thc first copy. Due to that, it is also not pos- 
sible to consider a purification protocol, where more than 
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just two copies are considered, and only CNOT-gates are 
required. 

Let us now show that the protocol Va, i.e. the purifi- 
cation of color A, achieves the following mapping: 

|*k A ,k A ) |*u,u) S kaM |* kA , kA ®i A ) . (3) 

Let us first investigate the projection onto the GHZ- 
states, li 2 © Wi \GHZ), with W t e {A^}^ k=0A . Dc- 

b 1 6 2 6 3 

noting by Pq H ' z z the projector onto the GHZ-st&te, we 
have l^^f = \ \GHZ) m (|0) bf <00| fc}fe? + 

|l>6f <H| 6 i 6? ) = \\GHZ) bim P b >. Note that 
P bi Wbl \4>)v> = O |0»63'. where |V>) © |0> denotes 
the Hadamard product between the states 1^) and \4>). 
That is |V> O 10) = Ei«iAI*> fOT W) = Ei«il*> and 
|0) = Ei A K)- As can be easily shown the projection 
onto the state 1 12 © VF, |GffZ) leads to VF, T P b >. The lo- 
cal unitary operation is undone by applying W* (see step 
(ii)). Note that the probability of measuring one of the 
four states {1 12 © W \GHZ) , with W G {X^Z 11 }^ k =o,i}, 
i.e. the probability of successfully implementing the pro- 
jector P bi is larger or equal to 1/2, as can be easily seen 
as follows. Let pi — J2 nm r nm \n) (m\ denote the reduced 
state of particle hi of p. Then, the success probability, 
Psucc = tr(P fc *p l © Pi P bi ) = tr{ Pi Q Pi ) = rl + r\ x . Us- 
ing that tr(pj) = r 00 + r n = 1 we find pf ait = l-p SUC c = 
2r 00 rii. Since p succ - p fa u = (r 00 - r n ) 2 > we have 

Psucc ^ 1/^- 

Let us in the following denote by P Na * the prod- 
uct of P b i for all neighbors, bj, of a^. Note that 
P Na z \tp) Uph |+)®' A ' = U p h for arbitrary phase gates 
Uph and |A|-qubit states \ip). 

With all that it is now straightforward to prove Eq (3). 
Combining steps (i) and (ii) we have: 

|$) - ® P N °> (g) ug'° l \^ aMa ) |* 1a , 1a ) = (4) 
(g) P»>< Y,*a\ \*^Ma) Ai , Ai Ha) (U\ A2 Z% 

ai£A j A 

m A 

Using now that both, <S> a . eA P Na * and U mA only act on 
the neighbors of A\, A 2 resp. and that Z£ j^k^.k^) = 
|*k A ,k A ei A ) we have 

m = E(- 1 ) U j " X A 1 ® l*k A ,^®i^) |j^) A2 = (5) 

JA 

J2\^^i a )\u)a 2 (-1)^^ = 

|*k A ,k A ei A )P 8|A| |Uek,i). 

Note that the first equality in Eq (5) is due to the fact 
that X A A x © J7j A = riiG{i |A|} ) which implies that 



X A \®U iA |* kA ,k A ei A ) = (-l) k ^ |* kA ,k A ei A ). Mea- 
suring now the qubits af e A 2 in the X -basis and keeping 
only those instances where all measurement outcomes are 
+1 leads to Eq (3). 

The intuition behind the subprotocols is the following. 
The projection on the GHZ-st&te applied to the neigh- 
bors of qubit Oj, combined with the application of the 
CNOT gate on the qubits a},af for all i maps the in- 
formation about the parity \s.a © 1a into the qubits in 
A 2 , which are then measured. That is, some information 
about the first copy is mapped to the second, which is 
revealed by the measurement (see Eq (3)). This infor- 
mation gain leads to the purification of the state. As 
mentioned above the subprotocol Va maps an arbitrary 
input state p = J2k k' A k l^k) (^Vl to the unnormalizcd 
state p = £ k k , A k ' |* k ) (* k ,|, where 

\k' A M' A _ \ ■> \ k ' A ,k A . k' A ,l' A ©k A , . 

k A ,k A - A k A ,k A A k A ,l A ek A - \°) 

1aA' a 

In the particular case of a LMES-diagonal state, p = 
E k Ak l^k) (^k|, the output is the unnormalized LMES- 
diagonal state, p — J2k A k l^k) (^k| with 

Ak A ,k A = E A k A ,1a A k A ,i A ek A ■ (7) 

Note that this equation resembles the one obtain for the 
purification protocols for graph states even though the 
protocol itself is fundamental different. The reason for 
that is that in both cases the second copy of the state 
is used to obtain the information about the bit values 
© 1,4 corresponding to parties in set A and in both 
cases only those instances where k^ © 1^ = are kept. 

In order to illustrate the subprotocol we consider the 
simple example of a n = 3m-qubit 3— colorable regu- 
lar LMES |*) = n?^ 1 1+)®", for m being an 
arbitrary integer. In Figure 1 we depict the subproto- 
col Va for this state. The stabilizers of color A have 
the following form {5" ai = X ai © U blCl ,S a2 = X a2 © 

Ub 1 c-iU Cl b 2 Ub 2C21 S a3 — X a3 © Ub 2C2 U C2 b 3 Ub 3C3 7 . . . S am — 

X am © Ub m _ lCm _ 1 U Cm _ 1 b m Ub m c m } ■ As explained above, 
the information about k^ is transferred to the second 
copy by applying CNOT gates between the qubits in A\ 
and A 2 and by performing the GHZ-measurements on all 
neighboring qubits, bi,Ci, for i = 1, . . . m (see Figure 1). 
The information is then gained by measuring the qubits 
af E A 2 . 

B. Total Purification Protocol for regular LMESs 

The total multipartite entanglement purification pro- 
tocol is composed of the subprotocols, Vx, for all color 
X. First Va is applied to two copies of the initial state. 
Then, Vb is applied to two copies of the output state, 
etc. Note that the convergence of the protocol depends 
on the order of the applied subprotocols. For certain 
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Figure 1: This graphic illustrates the protocol for a 
3— colorable LME state where color A is purified, see Sec. Ill A. 
The target state is given by |* ) = IPS ^-l.i.i+i 
The parties in color A entangle copy 1 with copy 2 by apply- 
ing CNOT-gates. All other parties make a GHZ-measurement 
on one qubit of each copy and one qubit of a maximally entan- 
gled state, |$ + )- This measurement is depicted as a triangle. 
Then a measurement in the X— basis is performed on the 
qubits of color A of the 2nd copy. This measurement is dis- 
played by an arrow pointing in the upper right direction. The 
qubits of the remaining state are colored in a darker fashion 
(online: bright red). 



input states an adaptive scheme is more advantageous 
than an alternating application of the subprotocols, as 
we show in Sec.VA. 

Due to Eq. (6) the new coefficients, A k , obtained af- 
ter applying one subprotocol, are non-linear functions of 
the coefficients A k . Thus, it is in general hard to find 
conditions under which the purification protocol, which 
is composed of several subprotocols Vx , converges to the 
desired state. Nevertheless, it is easy to see that a neces- 
sary condition for convergence is that the fidelity with the 
target state, An, is sufficiently large. In Sec.V we deter- 
mine the threshold for the fidelity for which the protocol 
still succeeds for several different states and noise mod- 
els. However, one must keep in mind that, due to the 
complexity of multipartite entanglement, this single fig- 
ure of merit cannot be the only relevant quantity (e. g. 
one can find pairs of states with the same fidelity with 
the target state but one can be purified while the other 
cannot). 



C. Purification Protocol for 7r— LMESs 

In this section we explain how arbitrary 7r-LMESs 
can be purified. Let us call a phase gate, U p h, regu- 
lar phase gate if the corresponding LMES, U p h |+}®" 
is a regular LMES. First of all, we note that an ar- 
bitrary phase gate U p u can be decomposed into regu- 
lar phase gates, where the number of qubits the regu- 
lar phase gates are acting on can vary. For instance 
the gate U p h — C/123^23 is decomposed into the reg- 
ular phase gates, Ux2$ and t/23. Second, note that 
two LMESs, |*0 = (U$ h ) iu ..., ik and |* a ) = 

(U$h)h,.~*i \+f\ for l,k € V with {h,...,i k } C 
{ji, . . . , ji} can be transformed via LOCC into the LMES 



l*> = (^k,...A(^)ii,..J. l+> This can be easily 
achieved by applying the GHZ-measurement explained 
in step (ii) in Sec. Ill A among all qubits in . . . ,if.}. 
Note that the condition {ii, . . . , i^} C . . . , ji} is not 
a restriction since one can always add qubits to the set 
{jx, . . . ,ji} onto which U^ h is acting trivially. Thus, in 
order to purify an arbitrary LMES, one would first de- 
compose the phase gate into regular phase gates, purify 
each of the corresponding regular LMESs and then com- 
bine the regular LMESs to obtain the desired state. For 
certain states there are definitely more suitable proto- 
cols than the one outlined above. For instance, the state 
p = x\V){^\ + i^l, with I*) = C/123C/23I+)® 3 can 
be purified to |\&) as follows. First, measure the first 
qubit in the Z-basis. In case outcome +1 is obtained, 
i.e. the state of the first qubit is projected onto |0), the 
state is (up to normalization) p = x {^2\ + ^pl) 
where 1*2) = ^23 I+)* 82 - This state can then be pu- 
rified to 1*2), which can then be used to perform the 
GHZ-measurement between p and 1*2) (^2! to obtain 

pt = pO(I*2) (* 2 |) = x |*i) (i-il + ±=21, with = 

^123 |+)® 3 - Since |* x ) is a regular LMES p' can be puri- 
fied using the methods explained in the previous sections. 
In a last step the GHZ measurements on one copy of |*i) 
and one copy of 1*2) are performed to obtain |\T/). 

IV. DEPOLARIZATION 

In this section we study how one could bring an arbi- 
trary state p — ^ k j A k |*k) ('I'll to its diagonal form in 
the LME basis p d = £) k A k |* k ) (* k |. For LME states 
this cannot be achieved in general by a deterministic 
LOCC map as we will show in this section. This is in 
contrast to the stabilizer states and graph states where 
the process can be implemented locally and where the 
depolarization is used as a first step in purification pro- 
tocols since this simplifies the analysis. We will never- 
theless show later how this can be done non-locally using 
LME states as a resource. 

In order to depolarize an arbitrary state p to an LMES- 
diagonal form, one can apply the map 

Stwirl(p) = £n • • -£i(p) (8) 

with Si = i(p + SipSj) and the Si are the stabilizers of 
the LME state \^o). The map makes the off-diagonal 
terms vanish, since S z |* k ) (*i| Sj = (-l) fc -®^ |* k ) (*,| 
and thus, for hi ^ In the off-diagonal terms of p and 
SipSj cancel each other. Notice that -unlike in the case 
of graph states- the operations Si are non-local. 

In the following we are going to show that the depo- 
larizing map and physically reasonable variations of it 
cannot be implemented by deterministic LOCC already 
in the simplest case of the 3-qubit regular LME state 
of order 3, i. e. |*) = C/i23|+)® 3 . This raises a second 
question: depolarizing simplifies the analysis of the pu- 
rification protocol but it is not necessary as the scheme 



7 



can nevertheless be applied to any state. Therefore, is 
it always advisable to depolarize before applying the pu- 
rification map (in the sense that this enlarges the set of 
states whose fidelity with the target state can be brought 
arbitrarily close to 1 by the protocol)? In the bipartite 
case it is known that depolarization can actually spoil pu- 
rification, in the sense that distillable states are mapped 
to non-distillable states. Here we show that this is also 
the case for our multipartite protocol, where we provide 
a state p for which our protocol succeeds in purifying to 
the 3-qubit regular LME state of order 3 but fails for 
£-twiri{p)- Notice that, since the mathematical form of 
the purification map is the same, this conclusion could 
be extended as well to 2-colorable graph states. In sum- 
mary, the purification procedure to LME states must be 
applied to possibly undepolarized states since in this case 
they cannot in general be made diagonal in the corre- 
sponding LME basis by LOCC. However, this needs not 
be regarded as a drawback as depolarizing does not neces- 
sarily contribute to the success of the purification process 
and there are states for which depolarization would actu- 
ally spoil the protocol. We will conclude this section by 
discussing different ways to implement the depolarizing 
map (8) at the expense of consuming entanglement. 

A. Are depolarizing maps in the LME basis locally 
implementable? 

In this subsection we consider the use of a depolarizing 
map Edep to preprocess the noisy state before the appli- 
cation of the purification protocol. One possible way of 
depolarization using a probabilistic application of stabi- 
lizer operators was already outlined above, however due 
to the non-locality of the stabilizers the corresponding 
operations cannot be implemented by LOCC. This does, 
however, not exclude the possibility that other sequences 
of LOCC can achieve the desired goal. 

The most reasonable assumption for £d ep is that for any 
input state the off-diagonal entries are suppressed and the 
diagonal entries remain equal. Moreover, the map should 
be independent of the state to be depolarized as we are 
aiming to devise general purification procedures which 
assume no previous information on the form of the noisy 
states to be purified. In other words, for any n-qubit 
state p it should hold 

(*k|£fep(p)|*k) = (*k|p|*k) Vk, 

(9 i \£ dep {p)\* k )=0 Vj^k. (9) 

Let 

A k = |*k)(*k| 

B Jk = |*j><* k | + |* k )<*j| (j<k) 

C jk = i(|*j><* k |-|* k )<*j|) (j<k), (10) 

where the relation j < k is to be understood to hold when 
j and k are mapped from binary to decimal basis (i. e. 



in lexicographical order). It is straightforward to check 
then that 

£ dep (A k ) = A k , £ dep (B ik ) = 0, £ dep (Cj k ) = Vj < k. 

(11) 

Since {A kl Bj k ,Cj k } is a basis of the real vector space 
of Hermitian matrices, the map £dep fulfilling Eqs. (9) is 
then unique for density matrices and therefore £d ep = 
£twiri- Now, if the corresponding Choi-Jamiolkowski 
(CJ) state associated to a map £ is entangled, then 
the corresponding operation cannot be implemented by 
LOCC [22]. Therefore, our question can be decided by 
computing the CJ state associated to £dep for different 
LME states following Eq. (8) or Eqs. (11). In the sim- 
plest case of the 3-qubit regular LME state of order 3 (i. 
c. S 1 =X 1 ® U 23 , S 2 = X 2 ® U13, S 3 = X 3 <g> U 12 ), one, 
more precisely, needs to compute the CJ state 

EiZ b2a2b3a3b = SSep ® l fc (l^ + }lal^ + )2a 26 |^ + ) 3a3fc ) 

(12) 

and check whether or not it is fully separable in the split- 
ting X 1 2 1 3. After lengthy but straightforward calculation 
one finds that 1 1 (-E123 ) r 1 1 tr — 1-75, where the superscript 
T indicates partial transposition with respect to one of 
the subsystems 1, 2 or 3 (the state is invariant under ex- 
change of these subsystems). Since ||i? r || tr = 1 iff the 
state has positive partial transposition, the correspond- 
ing CJ state i^lf is not fully separable and, therefore, 
£dep cannot be implemented by deterministic LOCC in 
this case. 

One may nevertheless wonder whether a wider class of 
maps £' dep with less stringent conditions but which still 
output states diagonal in the LME basis is still amenable 
by deterministic LOCC. Since one of the most relevant 
quantities for the purification to be possible is the fidelity 
of the noisy state with the target state, one may relax 
conditions (9) to 

<*olO)l*o) = <*oH*o>, 

(*jl^ep(p)l*k>=0 Vj^k, (13) 

for any n-qubit state p [26]. In addition to this, we fur- 
ther restrict £' dep to leave invariant all states which are 
already diagonal in the LME basis. We will show, how- 
ever, that the only positive map (i. e. mapping positive 
semidefinite matrices to positive semidefinite matrices) 
fulfilling the above conditions is again £d ep — £twirl and, 
hence, nothing is gained with such a relaxation. To see 
this, notice that the aforementioned restrictions on the 
action of £' dep boil down to 

£' dep (A k ) = A k Vk (14) 

and, together with the fact that £' dep should be trace- 
preserving, to 

^ep(Sjk) = E< A P ( 15 ) 
p#0 
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with J2 P = Vj,k and similarly for Cjk- Clearly 
then, for each j and k at least one of the real numbers 
a p k must be negative for some value p (if £' dep {B^) ^ 0). 
However, using linearity and for some fixed choice 1 we 
have that 

£' dep (aA + bAi + cBoi) = aA + bA x + c ^ a p l A p (16) 

for any a, b, c £ R (c > 0). Now, if a" 1 is negative for 
some p 7^ 1 the map is clearly not positive so af 1 must be 
negative. However, for any such possible negative value 
one can choose b and c such that b + aj" c < (i. e. 
the output of the map is not positive semidefinite) and 
still choose a such that ab — c 2 > (i. e. the input to 
the map is positive semidefinite). Thus, one sees that 
in this case the map is also not positive. Therefore, the 
only possibility is ^ ep (-Bjk) = Vj ^ k and obviously the 
same argument leads to ^ ep (Cjk) = Vj ^ k. Thus, this 
leads to £dep, which cannot be implemented by LOCC in 
general for LME states. 

For the sake of completeness, let us finally consider the 
most general class of maps that output states diagonal 
in the LME basis preserving the fidelity with the target 
state. That is, those maps $ satisfying just Eqs. (13) 
(dropping the condition that they leave invariant states 
which are already LMES-diagonal). We now have that 

<&dep(A ) = A 

t>dep(A k ) = PkmA m Vk ? (17) 

with P a (right) stochastic matrix (i. e. with all entries 
non-negative and ^ m Pkm = 1 Vk) due to trace preser- 
vation. For the same reason as above, one can see that 
^dep(Bju) = $dep(Cjk) = Vj ^ k must hold (since for 
any < P\ m < 1 and any a^J < one can choose b and 
c such that P\ m b + a^c < and still choose a such that 
ab — c 2 > 0). Hence, the support of & d e P is in the sub- 
space spanned by the {^4k} and without loss of generality 
it can be considered as a composition <& d ep = $p o £ dep 
where $p acts as given by Eqs. (17) for some choice of 
stochastic matrix P. This means that the CJ state for 
any $ dep , E 9 , is given by 

E? aU ... nanb = $ P <8> t b {E d Z b ... nanb ), (18) 

but this does not necessarily imply that all §dep cannot 
be implemented by LOCC since $p for some choice of 
P might break the entanglement of the CJ state of Sdep- 
However, for the 3-qubit regular LME state of order 3 
the P we have found decreasing the entanglement of E^ 
the most is still such that IK-E^^Htr = 1-165, which 
implies that Ef 23 is nevertheless entangled. Moreover, 
we have generated 1000 maps $dep with random choices 
for P and found that they could not be implemented by 
LOCC as the corresponding CJ state remains entangled 
(see Figure 2) . Therefore, even in the most general case of 



depolarizing maps that just preserve the fidelity with the 
target state, there is strong numerical evidence indicating 
that no such map can be implemented by LOCC for the 
simplest case of the 3-qubit regular LME state of order 
3. 

1.32 
1.3 
,_ 1.28 
U 1.26 



1.22 
1.2 
1.18 
1.16 



i 

Figure 2: |(-Ef23) r ||tr for randomly generated maps <fr' dep 

(1 < i < 1000). In all cases it holds that ||(£;f2 3 ) r ||tr > 1, 
indicating that the CJ state is entangled and, therefore, that 
the corresponding map cannot be implemented by LOCC. 



B. Depolarization may spoil purification 

As stated above, in many purification protocols a de- 
polarization to the basis corresponding to the target state 
is performed before the actual protocol is applied. This 
simplifies the mathematical analysis and one might think 
that this is favorable for the convergence of the purifica- 
tion process. By convergence we mean that the target 
state can be reached up to any given accuracy by ap- 
plying the protocol sufficiently often. However, in this 
subsection we construct a counterexample which shows 
that off-diagonal elements can improve the threshold of 
convergence. Similar results are known in the bipar- 
tite case. Consider, for instance, the case of purifying 
noisy states p = A^ n |*y fc )(* ;mn | to the 3-qubit 

3 -colorable state |* oo) = ^123 \+}® 3 - Using Eq. (6), one 
can identify conditions on the off-diagonal entries {A'™} 

such that they add to An,o but do not add to the value 
of the other diagonal elements no matter which color is 
being purified, thus obtaining a larger fidelity with the 
target state after one iteration of the protocol than if the 
state was depolarized. This happens to be the case if it 
holds that AjjJJ = -AgJ? e M, AjoJ = -A^? € R and/or 
Aqqq = — Aq^o S K and the rest of the off-diagonal entries 
are set to zero (up to the hermiticity condition Afe n = 




200 400 600 800 1000 
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(Kmn)*)- Therefore, one can find the desired counterex- 
ample by taking a diagonal state for which our purifica- 
tion protocol does not converge and populating strongly 
the off-diagonal elements as given above (as much as pos- 
itive semidcfinitcncss allows). For instance, consider the 

state p(f) with Ao,o = /, Ai,i = (I — f)/7 for i ^ and 
\on _ _\oio _ uoi _ _\ioo _ \iio _ _\ioo _ nni 

^000 ~~ A 001 ~~ ^000 — A 001 — A Q00 — A 01Q — u ' uz 

(with the other off-diagonal entries up to hermiticity be- 
ing zero). Notice that p(f) is positive semidefinite when 
0-01 < / < 0.72. Moreover, if / > 0.6503, the purifi- 
cation protocol converges. However, the twirled state 
£twiri{p{f)) can be purified only if / > 0.6507. Although 
the difference is small, this is enough to prove that it 
is not always more advantageous to depolarize the state 
before the purification protocol. 

C. Depolarization with nonlocal operations 

We have shown above that the map Stwiri given by Eq. 
(8) that depolarizes to an LMES-diagonal state does not 
belong in general to the class of LOCC operations. How- 
ever, this map can be implemented with certain proba- 
bility of success at the expense of consuming some entan- 
glement the parties might have initially at their disposal, 
or have generated from several noisy copies of their state 
with the help of entanglement purification protocols for 
different target states. One way to achieve this is the fol- 
lowing. For any map £i the parties decide randomly to 
apply Si or not. This leads to £i — Mp+SipSj). The im- 
plementation of Si = X % ®Ui can be easily done by apply- 
ing the GHZ-measurement that we have considered in our 
purification protocol (cf. Sec. Ill A), which succeeds with 
non-zero probability. In particular, the X of any Si can 
be applied locally by the ith party while the phase gate is 
implemented by using an entangled state \<j>) — Ui |+)® m 
shared among the m parties holding the qubits in Ni 
(to ease the notation we we do not explicitly write the 
dependence of |</>) on the LME basis). This is because 
Si |* k ) = Xi ® P Ni \<j>) |* k ), where, as before, P N * de- 
notes the product of GHZ projections on all neighbors of 
i. Graphically this is shown in Figure 3 a). Notice that 
Ui does not need to be pure and, hence, the m parties 
holding \<p) do not need to be fully entangled. Moreover, 
the phase gates which generate the auxiliary state | <j>) are 
of one order less than phase gates generating the LME 
basis in which the state is going to be diagonalized. For 
regular LME states generated by phase gates of order 3, 
Ui consists therefore of two qubit phase gates only and 
the state |</>) is then a (not necessarily connected) graph 
state for which purification protocols already exist. 

The map Stwiri can also be implemented by consum- 
ing pure LME states of the same order as the one 
that corresponds to the basis in which the state is go- 
ing to be depolarized. In order to implement the map 
Si one has to use an auxiliary state = (|0) (0^ ® 
1 + |1) ® Ui) |+)® ro+1 and apply the CNOT be- 
tween the qubit i, i' and P Ni on the neighbors to ob- 



tain £i(p) = tii,(U$'P Ni (p® {ip'\) P Ni Ug). Here the 
primed index i' is referring to the auxiliary state 
Figure 3 b) displays the procedure for one stabilizer Si. 




Figure 3: In this graphic we illustrate two methods to real- 
ize the CPM Si — ^(p + SipSj) of the depolarization process 
given by Eq.(8) where |*) = rL =2 , 4 ,6... ^-i.M+i 1+)®" Part 
a) depicts how Si is implemented using the auxiliary state \<j>) 
(see main text). The measurement in the GHZ basis is pic- 
tured as the triangle which acts on p, the auxiliary state, and 
one qubit of a maximally entangled state. The circle marks 
that a X operator is applied. In order to realize the CPM Si 
the parties choose randomly to apply this transformation or 
not. Part b) shows how Si is implemented by employing the 
LME state \ip'} and tracing out the qubit which is marked 
with a cross. 



V. EXAMPLES AND COMPARISON 

In this section we demonstrate the performance of our 
protocol by considering states, which are subjected to 
local as well as global noise. Then, we compare the di- 
rect purification protocol to indirect purification proto- 
cols, which, as we will show, can also be employed to 
purify to certain LMESs. 

A. Pauli channels 

As examples of the noise model we consider first the 
local depolarizing as well as the local dephasing channel. 
The CPMs describing the local channels are given by 
£{p) = ®j Si(p) where Si is either 

£i(p) = PP + ^-j^(p + X lP X t + YipYi + ZipZi) (19) 
for the depolarizing channel or 

£i(p) =pp + ^~^(p + z *P z i) ( 2 °) 

for the dephasing channel. In the following we will use 
£depo{p)(£deph(p)) whenever we refer to the depolarizing 
(dephasing) channel. We determined numerically the 
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smallest value of p for which the state still converges 
to the target state. Note that the convergence depends 
on the sequence of subprotocols we choose for the 
purification. Moreover, the most advantageous sequence 
depends on the input state. However, since we consider 
the case where one does not have any information about 
the noisy input state we investigate the examples for a 
fixed sequence. We consider up to 6— qubit 3— colorable 
LME states with colors A, B, C. We iteratively apply 
the order ABC — CAB — BCA. This sequence seems 
to perform better than just repeating the sequence 
ABC for a generic input state as the qubits in C get 
more noisy due to the backaction from the subpro- 
tocols purifying A and B. The states we investigate 
are given by U123 1+)® 3 and U123U234 \+) m for the 
3— and 4— qubit examples. For the 5 and 6 qubits 
we choose E/i23#234^345 \+)® 5 , ^123^124^125 I +) 85 

and ?7l23C/234C/345^ r 456 |+) 86 , ^134^235^234^346 I + )® 6 , 

[^123^124^125^126 I+)* 816 . The results we obtain for 
the depolarizing and dephasing channel (see Figure 4) 
suggest that there are three crucial properties which 
influence the tolerated noise level, namely the inter- 
action pattern of the target state, the colorability of 
the target state and the number of qubits of the state. 
The numerical results for the case of 3,4,5, and 6 
qubits for different states are shown in Figure 4. For 
the depolarizing channel states with a linear pattern, 
YYi=2 I+) 8 "*! seem to be more robust than the 

ones with GHZ— like interactions, e.g. Yl^ =3 U\2k \+)® n - 
The noise level does not strongly depend on the number 
of qubits. Note that similar results have been obtained 
for graph states [9]. Note further, that the threshold 
for p considering the depolarizing channel is close to 
the threshold of the dephasing channel even though 
the depolarizing channel outputs states which are non- 
diagonal and the dephasing channel outputs diagonal 
ones. This supports our previous observation, that it is 
not drastically advantageous for the purification process 
to have LMES-diagonal states. 

Let us now discuss how the colorability influences the 
noise thresholds. To this end we compare the tolerated 
noise level for the depolarizing channel acting on the con- 
sidered 3— colorable LMESs with the one obtained for 
the 2-colorable states, i.e. linear cluster states and GHZ 
states in [9]. We find that the 3— colorable states are 
more sensitive to noise. Note that this result is not very 
surprising, since, as mentioned above, the purification 
of one color adds additional noise to the other colors. 
Thus, if there are many colors in the state we need more 
rounds until one color is purified for the second time. 
Hence, it is more often subjected to noise before it is pu- 
rified again. From a physical point of view one might 
argue that a state with higher colorability is more con- 
nected and therefore the errors introduced by the chan- 
nels spread more easily. 

As a last example we also consider the influence of 
global noise. More precisely, we consider global white 
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Figure 4: In this figure we show the noise threshold p of 
the depolarizing and dephasing channel given in Eq.(19) and 
Eq.(20) and the minimal required fidelity for the purification 
process for a state subject to global white noise of Eq.(21). 
We consider 3— colorable LME states of 3, 4, 5 and 6 qubits 
which are generated by 3— qubit phase gates. The interac- 
tion pattern were given by Uviz |+)® 3 , C/123CV234 |+}® 4 for 
3 and 4 qubits respectively (x). For the 5— qubit state we 
have two possible patterns, a linear one U123U234U345 |+}® 
(x) and a GHZ-like one 

UimUimUvus l+>® 5 (o). In case 
of the 6— qubit state one has three different patterns, 
a linear one U123U234U345U456 |+}® 6 (x), an intermedi- 
ate one [/i34f235E/234f/346 |+)® 6 (o) and a GHZ-like one 
Ui23Ui2iUi2$Ui26 |+)® 6 (o). For the dephasing channel (□) 
the pattern does not strongly influence the threshold. For the 
global white noise channel (*) we show the fidelity for a target 
LME state with a linear pattern. 



denoting the number of qubits (see also Eq. (21)) acting 
on 3-colorable LMESs, |*) = [/;_i, M+ i \+)® n . We 

calculate numerically the threshold for the fidelity, / = 
($>\£ (I*) (*|) |*) such that the protocol still converges 
to the target state |\f r ). The numerical result for 3,4,5 
and 6 qubit states is shown in Figure 4, from which we 
can see how the threshold fidelity scales with the number 
of qubits in the state. 



B. Comparison to other purification protocols 

One way to generate a desired multipartite state is to 
generate maximally entangled bi-partite states shared 
between one party and all the others and to use tele- 
portation to distribute the multipartite state. Another 
approach would be to purify states, which can be used to 
generate (via LOCC) the desired LMES. For instance, we 
will show that it is possible to purify to certain LMESs 
by first purifying certain graph states and then recom- 
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bining them in order to obtain the desired LMES. If the 
performance of such a scheme would be better than the 
multipartite entanglement purification scheme, the later 
would loose its practical relevance. The aim of this sec- 
tion is therefore to compare the LME-protocol with pre- 
viously known purification protocols. We will consider 
some example to show that the LME-purification proto- 
col outperforms the graph state purification protocol as 
well as the bipartite purification protocol. 

Let us consider the 6— qubits linear LMES |\P) 6 = 

£^123^234^345^456 I+)* 816 - This state has the property 
that by measuring some of the qubits (of several copies), 
one can obtain certain graph (bi-partite) states, such 
that those states can be recombined via LOCC to again 
obtain the LMES. Thus, considering certain noise, one 
possibility to purify to the LMES would be to first ap- 
ply the measurements to obtain the graph (bi-partite) 
states; purify them via the previously known entangle- 
ment purification protocols and then recombine them to 
the LMES. In order to show how these schemes com- 
pare to the direct LME-purification protocol we consider 
states which correspond to the output of a global white 
noise channel, 



£ wn m (*|)=a;|*) <*| 



1 - x 



(21) 



where |*) = |*) 6 and x = f - ^r^r with / = 
(*|)|*) denoting the fidelity of the output 
state and n = 6 is the number of qubits. 

Before we explain several different schemes to purify 
to \*f>) indirectly, we state the main result of this sec- 
tion. We consider the situation where an arbitrary num- 
ber of copies of the state given in Eq. (21) is available 
and determine numerically the fidelity for which the var- 
ious purification protocols succeed. Applying the LME- 
purification protocol presented here, we find that the pro- 
tocol converges to \&) as long as the initial fidelity, /, 
fulfills / > 0.345. However, purifying the graph states 
which are needed to generate the LME state would re- 
quire a fidelity of / > 0.349. Using a bipartite strategy 
one requires a fidelity of / > 0.5 [5]. Therefore, using the 
protocol we introduced here, one can purify states which 
are too noisy to be purified with any other strategy. 

Let us now explain how the indirect purification pro- 
tocols work. We consider first the scenario where noisy 
graph states are obtained from Eq. (21) by performing 
local ^-measurements. Those states can then be recom- 
bined via LOCC to obtain the desired LME. We consume 
four copies of the state given in Eq. (21) to achieve this 
task. Let us denote by the measurement outcome 
of Zj measured on the i-th copy. First of all, the 2nd 
and 5th qubit of copy 1 and 3 are measured in the Z- 
basis. We keep only those instances where the outcomes 
m\ = ljTOg = — 1, and = —\,m\ — 1 are obtained. 
In the following we do not consider the qubits which 
have been measured, since they factor from the result- 
ing noisy graph (bi-partite) state. The resulting state of 



copy 1 (3) is given by Eq. (21) with |*) = Ui 3 U 34 
= U 3i U i6 \+f 3 ) resp. (see Figure 5 a). 
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a) 



b) 




Figure 5: This figure illustrates how one can obtain graph (bi- 
partite) states from LME states by local Z— measurements. 
The input state is given by £(|*} (*|) from Eq.(21) with 
I*) 6 = ^123^234^/3451/456 |+)® 6 - The local measurements in 
the Z basis are pictured by arrows pointing in the upper right 
direction and the numbers +1,-1 denote the outcomes. Part 
a) of the figure shows how one can obtain a 3— qubit graph 
state <*|) with j*) 3 = U13U34 \+)® 3 ■ Part b) indicates 

how a bi-partite state £(|*) (*|) with |*} 2 = U 12 \+) m is 
obtained. The coloring of the qubits refer to the three colors 
of the LMES. 



On the second copy qubit 3 and 5 are measured. The 



copy is kept only if m| = 1 and = — 1. The resulting 



®3 



state is then given by Eq. (21) with = U12U24 
Measuring qubit 2 and 4 of the fourth copy and obtain- 
ing m| = — 1 , m\ = 1 leads to the state given in Eq. 
(21) with |\&) = L/35C/56 |+)® 3 . After these measurements 
all graph states are purified using the purification proto- 
col presented in [8]. Note that all states are 2-colorable 
graph states. Note that the output of this process has 
the same fidelity with the target graph state as the orig- 
inal noisy LME state £ wn (\^) (*|) with the target LME 
state. Given the results presented in [8] it is easy to see 
that the protocol converges only if the original fidelity 
was / > 0.349. Next, we explain how these four pure 
graph states can be mapped to the desired LMES via 
LOCC. To this end let us denote by Qi the measurement 
operator pu>H ® H, acting on qubit i of two different 
states. Here, p i)V = V2(|0) (++| + |1) (++| U iti >, where 
Ui t i' denotes a 2— qubit it— phase gate acting on the i-th 
qubit of both states and H denotes the Hadamard op- 
eration. Whenever it is clear from the context we will 
omit the subscript. First, the states which were obtained 
from copy 1 and 2 are recombined by performing the lo- 
cal measurements Q\ and Q4 (see Figure 6 a). On copy 
3 and 4 Qi for i = 3, 6 is performed. Finally, the opera- 
tions P 3 and P 4 with Pi = v / 2(|0) i (00| ii , + |l) i (1%) are 
performed on the two LME states which were obtained 
in the previous step (see Figure 6c)). The resulting state 
is the desired LMES, |*). 

An other approach is to distill pure bi-partite maxi- 
mally entangled states and to use them to generate the 
LMES via LOCC. In order to get a noisy two qubit state 
one could measure 4 consecutive qubits in the Z-basis 
obtaining the outcomes +1,-1,— 1,-1 as illustrated in 
Figure 5 b). In order to purify maximally entangled 
two-qubit states an initial fidelity of / > 0.5 is re- 
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quired [5]. Once the bipartite states are purified one 
can recombine them again to obtain the LMES in the 
following way. We use two of the maximally entan- 
gled states, say |5 >+ ) 12 and | ) 2 , 3 and perform the lo- 
cal projector, Q2H® 2 to obtain the 3— qubit LME state 
Q2H® 2 |$+) 12 |$+) 2 , 3 = Ui 2S \+f 3 (see Figure 6). |*) 6 
is then obtained by taking four of those 3— qubit LME 
states and connecting them in a similar way. For in- 
stance, to obtain UijkUjw |+) 86 for some i,j,k,l' G V 
we perform the measurements on qubits j,j' and on 
qubits k,k', i.e. PjPkU ijk U fkn < |+)® = U tjk U jM > |+)® 6 
(see Sec. Ill and Figure 6). 



a) b) 




Figure 6: This figure shows how one could connect two graph 
states (bi-partite) states to obtain LME states [a, (&)] and 
how one can recombine LMES (c, d) . In part a) the LMES 
EA23C/234 |+)® 4 is obtained by applying the projector Qi (see 
main text) between graph states of the form U12U24, |+}® 3 and 
UvyUyA' |+}® 3 . In part b) QiH® 2 is applied to maximally 
entangled states resulting in C/123 |+)® 3 . Part c) and d) 

show how LMES can be recombined by employing the pro- 
jector Pi (see main text). The coloring of the qubits refers in 
all parts to the colors of the quantum state. 



VI. CONCLUSION 

In this paper we have devised a multipartite entangle- 
ment purification protocol for LME states with 7r-phase 
interactions. These constitute a large class of states with 
interesting physical and mathematical properties, which 
includes and generalizes well-known subclasses such as 
stabilizer states and graph states. In particular, we have 
proposed a purification procedure for the subclass of reg- 
ular LME states which, in turn, can be used as a build- 
ing block for protocols in which the target state is not 
regular. The main idea behind the protocol traces back 
to the Bell-state recurrence protocol: two copies of the 
noisy state are used to obtain some information about 
the eigenvalues of the stabilizers which is revealed by 
measuring (i. e. destroying) just one copy. However, 
LME states have very different mathematical properties 



(e. g. the stabilizers are not local) and have more com- 
plex interaction patterns (which affects the way in which 
quantum correlations are spread through the state) and 
it has been necessary to develop significatively different 
tools and techniques to achieve this task. Another differ- 
ence with previous purification protocols is that, as we 
have proven, depolarizing to an LMES-diagonal form is 
in general not possible by deterministic LOCC. However, 
we have shown that depolarization is not necessarily ad- 
vantageous for the success of purification protocols. We 
have also considered different examples to illustrate the 
performance of our protocol, determining the maximum 
amount of local noise under depolarizing and dephasing 
channels and the minimal fidelity with the target state 
under global white noise that allows for the convergence 
of the purification protocol. Remarkably, we have addi- 
tionally shown that our protocol purifying directly to the 
LME state has a better tolerance than indirect strategies 
that purify to several different graph states from which 
the LME state can be obtained locally (i. e. our proto- 
col converges for states which are too noisy for the other 
strategies) . 

It is worth pointing out that by convergence of our 
protocol we mean that the fidelity of the noisy state with 
the target state can be brought arbitrarily close to 1 by 
repeating our scheme sufficiently often. However, the 
fidelity only reaches 1 in the limit of infinitely many it- 
erations. This means that, strictly speaking, the yield 
Y of our protocol, which is defined as the ratio m/n 
between the number m of perfect copies of the target 
state obtained from n copies of the noisy state when 
n goes to infinity, is zero. This is also true for pre- 
viously proposed purification protocols for graph states 
[9, 10] and is not very relevant from the practical point 
of view. However, from the theoretical point of view it is 
interesting to note that these protocols can be modified 
using the hashing and breeding schemes [5, 23] to ob- 
tain a non-zero yield for states with very large fidelity 
with the target state. The idea behind these proto- 
cols is to regard the n copies of the depolarized noisy 
state P ® n = Ek(i>-k(™> A k<D • • • A k<Il ) |* k w)<*kW I®- • -® 
l^'k("))(^ , k(™) I as a statistical mixture. One then needs 
to identify the bit string k = (kW , . . . , k( n )) that was 
actually produced by the source taking into account that 
the probability of emitting each of the bits in k is en- 
coded in Ak- Once k is determined the state is identified 
to be l^kt 1 )) <8> • • • <8> l^kt"))) which can be transformed 
by local unitaries to |\&o)®"- Let us consider, for in- 
stance, the first bit of k, k ai , which is with probability 
p(k ai = 0) = Sk\fe ai ^o*a 3 -k ai ,kj and 1 witn probabil- 
ity p(k ai = 1) = J2k\k ai *ifca 2 -*aj ,kj- Since tne proto- 
cols allow to read the parity of the bits, in the asymp- 
totic limit of many copies one just needs to measure (i. 
e. waste) nS(p ai ) copies of the state to determine the 
first substring of k, k ai = (fci^, . . . , k^). Here, S(p ai ) 
is the Shannon entropy of the probability distribution 
{p(k ai = 0),p(k ai — 1)}. Since the parity checks of 
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bits in the same color can be done jointly this protocol 
yields Y = 1 — maxj S{p ai ) — max^ S{pf )i ) — ••• , which 
can be strictly larger than zero. Notice that for the pu- 
rification of LME states one can also read jointly the 
parity of bits in the same color. However, in our protocol 
this is achieved probabilistically by doing several GHZ- 
projections. In general, the probability of implementing 
each of these projections successfully is strictly less than 
one. Therefore, a correct reading of the parity of some 
substring of bits requires exponentially many copies, i. 
e. scales exponentially with n, thus rendering impossi- 



ble a modification as above of our protocol to obtain a 
non-zero yield. 
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